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Results are presented from a search for the decays D° — > K~tt + and D° — > K + ty~ in a sample of 
3.8 x 10 6 central Pb-Pb events collected with a beam energy of 158,4 GeV by NA49 at the CERN 
SPS. No signal is observed. An upper limit on D° production is derived and compared to predictions 
from several models. 

PACS numbers: 25.75.Dw 



The measurement of open charm production in heavy- 
ion interactions is of considerable interest because charm, 
due to its large mass, is predominantly created at the 
early stage of the collision when the energy density is 
large. Because of the hard scale involved, perturbative 
QCD (pQCD) calculations can serve as a baseline for the 
study of the production mechanisms and the dynamical 
evolution of charm in these collisions. 

At present, no direct measurement exists of open 
charm production in heavy-ion interactions at SPS en- 
ergies. The NA38/50 experiment has, however, observed 
a significant enhancement of di-muon production in the 
intermediate mass range of 1.5-2.5 GeV, compared to di- 



muon yields expected from the Drell-Yan continuum and 
semi-leptonic charm decays |l||. The origin of this en- 
hancement is presently not clear but can be explained by 
assuming that open charm production in central Pb-Pb 
collisions is about a factor of 3.5 times larger than pre- 
dicted by pQCD. This enhancement is currently under 
investigation by the NA60 experiment . 

A variety of models give very different estimates for the 
open charm yields at the SPS. For instance, for central 
Pb-Pb interactions at 158A GeV beam energy, a pQCD 
calculation based on Pythia predicts a yield per event of 
7V(D° + D ) = 0.21 (the centrality is here characterized 
by the number of participant nucleons N paTt = 400) . 
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In a yield of 0.5-0.6 cc quark pairs is calculated, based 
on the J/\I/ yield measured by NA50 and the statistical 
coalescence model (iV par t = 360). This translates into 
7V(D° + D°) m 0.4 if one assumes that about one third 
of the cc hadronize into D° and D , like in p-p interac- 
tions |5j. The ALCOR hadronization model |6j gives a 
much larger estimate of W(D° + D°) = 2.4 (A part « 350). 

An even larger yield of N(D° + D ) ss 6 is predicted by 
the statistical model of the early stage (SMES) which 
assumes charm equilibration in a deconfined quark-gluon 
plasma (QGP) at the early stage of the Pb-Pb interaction 
(A^t =360) 0- 

To discriminate between the different model predic- 
tions and to possibly identify the origin of the di-muon 
enhancement seen by NA38/50, we have performed a 
search for open charm, using invariant mass reconstruc- 
tion, in a large data sample of about four million central 
Pb-Pb events collected at 158A GeV beam energy. 

The NA49 detector [8j is a large acceptance fixed- 
target hadron spectrometer at the CERN SPS. Track- 
ing is performed by four large-volume TPCs. Two of 
these are placed one behind the other inside two super- 
conducting dipole magnets (vertex TPCs). The two 
other (main) TPCs are placed downstream of the mag- 
nets left and right of the beam line. These main TPCs 
increase the lever arm of the track reconstruction and are 
optimized for particle identification through a measure- 
ment of the specific energy loss (dE/dx) with a relative 
resolution of about 4%. The combined TPCs provide an 
accurate measurement of the particle momenta with a 
resolution of Ap/p 2 w 3 x 10~ 5 (GeV/c) -1 . Centrality 
selection is based on a measurement of the energy de- 
posited by the projectile spectator nucleons in a forward 
calorimeter. 

To measure rare particles like the f2 [t| and to search 
for open charm a large sample of central Pb-Pb events 
was taken in the year 2000 with a beam energy of 
158A GeV. In this run 3 x 10 6 events were collected with a 
centrality selection of 23.5% of the inelastic cross-section 
(■Wpart = 262). Also included in the present analysis is 
a 1996 data set of 8 x 10 5 Pb-Pb events, taken at the 
same beam energy but with a 10% centrality selection 
(A/part = 335). To increase the data acquisition speed 
and decrease the data volume only every second time- 
sample of the TPCs was read out during the 2000 run 
(256 instead of 512 time-samples). The reduced sam- 
pling did not significantly affect the track reconstruction 
nor the dE/dx measurement. 

The D° were reconstructed via their charged parti- 
cle decays D° — > K~tt + and D° — > K + 7r~ (4% branch- 
ing ratio). Because the secondary vertex resolution of 
about 1 cm is not sufficient to detect the decay vertex 
(7CT rs 1 mm), the D° candidates were identified by 
selecting a window around the D° mass in the invari- 
ant mass spectrum of the daughter particles. With a 
multiplicity of approximately 1400 reconstructed charged 
tracks, about 5 x 10 5 entries for each event were made 
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FIG. 1: Invariant mass distributions of D° (top) and D can- 
didates (bottom) from track samples without (left) and with 
(right) kaon identification. The open (shaded) histograms are 
before (after) applying the decay angle cuts described in the 
text. The vertical lines indicate the D° mass window used in 
the analysis. The distributions are corrected for acceptance 
and branching ratio. 



in each of the D° and D invariant mass spectra lead- 
ing to a large combinatorial background. Because of the 
large multiplicities it was not possible to measure the 3- 
particle decay D* — > K7T7T even though the background 
in this channel is suppressed by kinematic constraints. 

Events for which the primary vertex could not be de- 
termined were discarded from the analysis. Several track 
quality cuts ^(| were applied to remove non-vertex or 
badly reconstructed tracks. The remaining sample was 
sub-divided into two classes. The tracks in the first class 
have sufficient track length in the main TPCs and low 
enough momentum so that a significant enrichment of 
the kaon content could be achieved by suitable cuts on 
dE/dx. For tracks in the second class this kaon identifi- 
cation was not possible. Loose cuts on dE/dx (2a around 
the Bethe-Bloch curve) were applied on the tracks in the 
first sample to minimize the loss of kaons (and D°). The 
identified kaon tracks were then combined with all oppo- 
sitely charged tracks and the invariant mass of the pair 
was calculated assuming that the associated track was a 
pion. In the second class (without dE/dx) the invariant 

mass oftheD (D ) candidate was calculated for all pairs 
of oppositely charged tracks assuming that the positive 
track was a pion (kaon) and the negative track a kaon 
(pion). The invariant mass distributions obtained from 
the D° and D samples with and without kaon identifica- 
tion (corrected for acceptance and efficiency, see below) 
are shown by the open histograms in Fig. ^ 
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FIG. 2: The reconstructed decay angle distributions (see 
text) of simulated D° decays (left) and of D° candidates from 
real events (right). In both plots the D° candidates have a 
reconstructed mass within 50 MeV of the nominal D° mass 
and a reconstructed pt of 800-1200 MeV. A darker shade cor- 
responds to more entries in the plots. The full lines indicate 
the cuts used in the analysis to discard the hatched regions. 



To further reduce the combinatorial background, decay 
angle cuts were applied as follows. For each D° candi- 
date the polar angle 9 and the azimuthal angle </> of the 
kaon track were calculated in the rest-frame of the D°. 
Here 9 is the angle between the beam direction and the 
kaon flight direction and <f> the angle between the kaon 
and the flight direction of the D° in the plane perpen- 
dicular to the beam. In the left-hand side plot of Fig. [21 
is shown the distribution of decay angles from simulated 
D° decays (see below). The distribution from real events 
(almost entirely background) is shown in the right-hand 
side plot. It is clear from this figure that the signal 
distribution is approximately flat while the background 
distribution peaks at large values of |cos0|. Cuts, like 
those shown in the figure, were optimized to maximize 
the significance (= signal / ^/background) of the measure- 
ment. Because the decay angle distribution depends on 
the transverse momentum (px) of the D° and is different 
for the samples with and without particle identification, 
separate cuts were determined, for each of the two sam- 
ples, in five pt bins of 400 MeV width. The decay angle 
cuts reduced the background by a factor of about 3 (10) 
in the sample with (without) particle identification. (The 
signal is reduced by 30-40%.) Cuts on other kinematic 
variables like rapidity (y) and px were investigated but 
were found to be ineffective in the separation of signal 
and background [Toj . The invariant mass distributions 
after the decay angle cuts and corrected for acceptance 
and efficiency (see below) are shown by the shaded his- 
tograms in Fig. ^ 

To determine acceptance, efficiency and mass resolu- 
tion a Monte Carlo sample of D° and D° mesons was gen- 
erated with a Gaussian distribution in y (a y = 0.6) and 
an exponential distribution in transverse mass (300 MeV 
inverse slope parameter). The D° (D ) and their decay 
particles were transported through the NA49 detector 
geometry using GEANT 3.21 [HI, followed by a detailed 
simulation of the TPC response using dedicated NA49 
software. The simulated raw data were added to real 
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FIG. 3: Invariant mass distribution of simulated D° de- 
cays embedded in real data events. The full curve shows the 
Cauchy parameterization described in the text. 

events and subjected to the same reconstruction pro- 
cedure as the experimental data. The acceptance was 
calculated in bins of y and px as the fraction of D° 
(D ) which are geometrically accepted, survive the re- 
construction procedure and pass the analysis cuts. The 
experimental acceptance covers the range px > and 
— 1 ^5 V 1-6 and is found to be, on average, 8.4 (12.0)% 
for the sample with (without) kaon identification. It 
was verified that the amount of accepted particles var- 
ied by only 10-20% if reasonable alternatives (e.g. from 
Pythia [lj]) were chosen for the kinematic distribution of 
the D°. The invariant mass distributions shown in Fig.Q] 
are divided by the acceptance and by the branching ratio 
for D — > 7rK decay. 

The simulated data served to determine the shape of 
the invariant mass distribution of reconstructed D° as 
shown in Fig. The shape can be well described by a 
Cauchy distribution (curve in Fig. |3) 

dn _ N T 

dm" ~~ 2^ (m- m Q ) 2 + (r/2) 2 ' ^ ' 

where N is the total D° yield per event, mo the D° mass 
and r the width of the distribution. This width is al- 
most entirely determined by the detector resolution and 
is found to be T = 6.2 MeV with the mass of the D° set 
to m = 1864.5 MeV 0. 

The invariant mass spectra shown in Fig. ^ were fitted 
(by X 2 minimization) in a region of ±90 MeV around 
the nominal D° mass to the sum of a signal distribution, 
Eq. QJ, and a fourth order polynomial describing the 
background. The position mo and width T of the signal 
distribution were kept fixed to the values given above 
while the normalization N was left a free parameter in 
the fit. This fit results in yields (per event) of N(D°) = 
-0.41 ± 0.51 and N(~D°) = 0.05 ± 0.54, where the errors 
are statistical only. In Fig.0]is shown the invariant mass 
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FIG. 4: Invariant mass distribution of the D° + D° candidates 
after background subtraction. The errors are statistical only. 
The full curve shows the fit to the signal distribution described 
in the text, the dashed curve the expectation from ALCOR 6] 
and the dotted curve that from SMES 0. 



distribution of the D° + D candidates after background 
subtraction. Clearly no signal is observed. The fit gave 
for the total yield a value of N(D° + D°) = -0.36 ± 0.74 
per event (full line in Fig. 0}. 

An upper limit for the number of D° per event is esti- 
mated in a Bayesian approach ^4|. Here the likelihood 
P(data|iV) (i.e. the conditional probability density dis- 
tribution of the data, given N D° per event) is parame- 
terized as a Gaussian 



P(data|iV) = 



1 



exp 



2a 2 



= g(N;n,a) 



with mean /i = —0.36 and width a = 0.74 as obtained 
from the \ 2 fit- Using Bayes' theorem the posterior dis- 
tribution P(iV|data) is calculated by multiplying the like- 
lihood with an assumed prior probability distribution of 
N which is taken here to be zero for N < and uniform 
for N > 0. This prior distribution forces N to be positive, 
as it should be. Integration of the posterior distribution 
gives for the confidence level 



CL = 



P(AT|data) dN 



f™g(N;»,<r) dN 
f™g(N;»,a) dN' 



(3) 



where M is the upper limit of N corresponding to the 
confidence level CL. The denominator on the right-hand 
side of Eq. JSJ) accounts for the proper normalization of 
P(iV|data). Using the fitted values of and a, the upper 

limit for the total yield is found to be M(D° + D°) = 1.5 
per event at 98% CL. 

Because no D° signal has been observed it is not pos- 
sible to directly verify the Monte Carlo prediction of the 
signal shape. To investigate the sensitivity of the upper 
limit to the width of the mass peak the fits were repeated 
with T = 12.4 MeV. This resulted in N(D°) = -0.46 ± 
0.85, N(D°) = -0.22 ± 0.90, N(D° + D°) = -0.7 ± 1.2 
and an upper limit of M(D° + D ) = 2.4 per event at 98% 
CL. We remark that increasing the Monte Carlo estimate 
of the width by a factor of two should be considered a 
very generous error on T. 

Taking as a standard for comparison the pQCD esti- 
mate of iV(D + D ) = 0.21 mentioned in the introduc- 
tion, we conclude that an enhancement of charm produc- 
tion by more than a factor of 5-10 at the SPS is very 
unlikely. Due to the large combinatorial background it is 
not possible to confirm, nor exclude, a charm enhance- 
ment by a factor of three allowed by the NA38/50 mea- 
surement. However, the D° upper limit from this analysis 
is only marginally compatible with the yield estimated by 
the ALCOR model (dashed curve in Fig. 0} and clearly 
incompatible with the equilibrium yield of charm in a 
QGP as predicted by the SMES (dotted curve in Fig.Q). 
The latter observation does not necessarily exclude QGP 
formation at SPS energies provided that the QGP life 
time is shorter than the equilibration time of charm. 
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